Non-extensive statistics in stringy space-time foam models and entangled meson states 
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The possibility of generation of non-extensive statistics, in the sense of Tsallis, due to space- 
time foam is discussed within the context of a particular kind of foam in string/brane-theory, the 
D-particle foam model . The latter involves point-like brane defects (D-particles) , which provide 
the topologically non-trivial foamy structures of space-time. A stochastic Langevin equation for the 
velocity recoil of D-particles can be derived from the pinched approximation for a sum over genera in 
the calculation of the partition function of a bosonic string in the presence of heavy D-particles. The 
string coupling in standard perturbation theory is related to the exponential of the expectation of 
the dilaton. Inclusion of fluctuations of the dilaton itself and uncertainties in the string background 
will then necessitate fluctuations in Qs. The fluctuation in the string coupling in the sum over genera 
typically leads to a generic structure of the Langevin equation where the coefficient of the noise term 
fluctuates owing to dependence on the string coupling gs- The positivity of ga leads naturally to 
a stochastic modelling of its distribution with a x~ distribution. This then rigorously implies a 
Tsallis type non-extensive or, more generally, a superstatistics distribution for the recoil velocity of 
D-particles. As a concrete and physically interesting application, we provide a rigorous estimate of 
an oj-like effect, pertinent to CPT violating modifications of the Einstein-Podolsky-Rosen correlators 
in entangled states of neutral Kaons. In the case of D-particle foam fluctuations, which respect the 
Lorentz symmetry of the vacuum on average, we find that the t<;-effect may be within the range of 
sensitivity of future meson factories. 



> 

a^ 
m 

00 

O 

> 

X 
J3 



I. INTRODUCTION 

The standard model of particle physics is considered to be successful even though it has many undetermined 
parameters which need to be fitted to data. One reason for the great interest in physics beyond the standard model 
(SM) is the absence, within its framework, of quantum gravity (QG). A full picture of QG still remains elusive. 
One interesting approach to it, based on string theory (ST) [l|, has the advantage of also unifying QG with 
the nuclear and electromagnetic forces. Moreover flavour mixing phenomena are not understood in the sense that 
phenomenological mixing matrices are not derived from a fundamental point of view. There may be some relationship 
between these two deficiencies, and QG may play a role in that. For instance, (a small) part of the observed mass 
differences between neutrinos might be quantum gravitational in origin as argued recently ^ within the context 
of a space-time foam model in string theory with space-time defects Q in the form of point-like D(irichlet)-branes 
(D-particles) The interaction of the defects with string matter in such a model induces fiavour mixing and 

non-trivial contributions to the cosmological constant or better vacuum energy) proportional to the mixing angle and 
the mass differences among the neutrino states. 

Because of its foundations based on local relativistic unitary quantum field theory, SM has CPT symmetry 0]. 
Hence the detection of possible violations of CPT necessitates clearly physics beyond the standard model. CPT 
violation might be an important feature of QG, which may be exhibited by some stringy models of space-time foam 
with defects, of the type discussed in Q. Recently, an interesting signature of CPT Violation due to decoherence of 
matter as a result of QG effects in stochastic space-time foam models, has been suggested in ref. 0, implying new 
physics beyond SM in the entanglement properties of neutral meson pairs. The proposed signature is the so called 
omega effect where modifications to the standard Einstein-Podolsky-Rosen (EPR) correlations appear for entangled 
states of neutral flavoured mesons in meson factories [S, 9,]. 

Neutral mesons, such as the K mesons, have in the past been pivotal in the study of discrete symmetries p^. 
The decay of a (generic) meson (such as. the (j) meson produced in collisions of and with quantum numbers 
J^'^ — 1 [Hi), leads to a pair state |i) of neutral mesons (M) which has the form of the entangled state 

|^) = (|llo ("fc ) \ mJ-^\)- Mn(~k]) Wj-~^ 



(1.1) 



This state has CP = If CPT is not well-defined (actually perturbatively, i.e. although the concept of the anti- 
particle still existsit is slightly modified), then Mq and Mq may not be identified and the requirement of CP = + 
can be relaxed @|,@. Consequently the state of the meson pair can be parametrised to have the form 
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where uj = |a;|e'^ is a complex CPT violating (CPTV) parameter [g. It turns out to be difficult to generate this 
Lo term since reasonable attempts based on local effective lagrangians with explicit CPT violating terms do not lead 
to it. Recent approaches have concentrated on properties of space-time foam and related decoherence, which by 
their nature go beyond the concept of local effective field theory. From extrapolations of black hole physics and 
quantum uncertainty relations to the smallest scales it seems unavoidable that space and time should fluctuate on 
a space scale Zp, the Planck length, and a time scale tp, the Planck time. Such quantmn fluctuations lead to a 
background state which will be generically called space-time foam .12]. The full understanding of quantum states 
of matter in such a medium foam requires a deeper understanding of QG. Consequently the approaches to it have 
ranged from the purely phenomenological to those based on a more fundamental theory such as strings. Even in the 
latter, phenomenological assumptions will have to be made since we do not have a developed string field theory and 
formulations have typically been based on fixed backgrounds. Indeed a plausible phenomenological ap pro ach based 
on thermal baths, at least in its natural and most straightforward form, cannot generate an w term [13|. On the 
other hand a stochastic approach based and inspired by theoretical considerations of D-particle recoil from stringy 
matter in D-particle foam models [3| has predicted and quantified this omega effect Q. However, the quantification 
was based on somewhat phenomenological and naive estimates. In this work we shall provide more rigorous estimates 
of this effect within D-particle foam models. 

In the treatments so far of D-particle foam the stochastic modelling has concentrated on the randomness of the 
recoil velocity vector of the D-particle defect during its topologically non-trivial interaction with the matter string 
state (i.e. the capture and subsequent re-emission of the matter state) . In this work we want to incorporate the effects 
of the (target-space) quantum fiuctuations of the D-particle which, within the bosonic string model, can be calculated 
using perturbation theory. As we shall discuss, the renormalization of a certain subleading divergence to all orders 
in string perturbation theory leads to fiuctuations superposed on a drift velocity of the D-particle obtained after the 
initial recoil on interacting with stringy matter. Furthermore, at the phenomenological level uncertainty in the string 
vacua can lead to a fiuctuating string coupling which is related to vacuum expectation value of the dilaton field. 
Incorporation of such fluctuations lead to superstatistics [3] of the velocity recoil distribution, non-extensive statistics 
of Tsallis being one possibility. The recoil fluctuations of the D-particle lead to effective stochastic back-reaction 
on space-time that cannot be neglected. Hence the space-time metric will have induced stochastic contributions from 
stochasticity in the recoil p^ . 

In this work we shall concentrate our discussion on the Bosonic string. Although admittedly, this is not relevant 
phenomenologically, it is the only case where we manage to sum over world sheet genera and express our basic 
results on stochastic fluctuations of space-time and fuzziness of string coupling in a closed form. We shall make some 
comments on the robustness of our results and thus their extension to supersymmetric cases at the end of our article. 
This constitutes an active part of our research at present. 

Within the Bosonic c-model framework of D-particle foam, we shall calculate the stochastic fluctuations of the 
space-time induced by D-particle recoil, and then we shall use it to estimate the resulting w-effect in the initial 
state of entangled mesons in this concrete microscopic model. Contrary to previous naive estimates @, based on 
dimensional analysis, in this particular example of (string-inspired) foam, the effect is calculated rather rigorously, 
following conformal fleld theory methods in the world-sheet of the string. In particular, for the magnitude of the 
w-effect in the initial entangled meson states (|1.2I) . we find the following result for the square of the amplitude of the 
parameter u: 

' ' (mi - TO2)2 Ml ^ ' 

for non relativistic entangled states of mesons, with (near degenerate) masses ,? = 1,2 and Mp the quantum 
gravity scale, assumed to be the four dimensional Planck scale, 10^^ GeV. For comparison, we remind the reader 
that the naive phenomenological estimates of lead to effects of order l^^^l naive ^ m'^ (ml-ma)^ ' ^^^^ ^ typical 

momentum transfer variable, during the interaction of the meson state with the D-particle, C ^ which in Q has 
been assumed not very small. What the detailed string calculation in this paper shows, in which the stochastic space- 
time fiuctuations are due to quantum fiuctuations of recoil- velocities about a zero (Lorentz invariant) average, is that 
the parameter C of depends on details of the meson system, such as the mass and momenta of the mass eigenstates, 
and is effectively of order ~ (rn\ -\- 'm'^)/k'^ . For the case of neutral kaons in 0-factories, with k ~ 0(1 GeV), this 
is already at the order of magnitude required for detectability in an upgrade in DA<I>NE [9,] . Hence these models of 
D-particle foam may be falsifiable in such future meson factories. Of course, the estimate assumes the neutral Kaon 
as interacting with the space-time foam as an elementary entity, ignoring details of strongly interacting constituents 
in the Kaon substructure. Taking into account such strong interaction effects might change the above estimates. 

The structure of the article is as follows: in the next section we will discuss general features of the D-particle foam 
model, on which our work is based. In section IIIII we review the basic mathematical properties of the deformations 
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FIG. 1: Upper picture: A fluctuating D-particle is described by open strings attached to it. As a result of conservation of string 
fluxes [2,[a,lB| that accompany the D-branes, an isolated D-particle cannot occur, but it has to be connected to a D-brane world 
through flux strings. Lower picture: World-sheet diagrams with annulus topologies, describing the fluctuations of D-particles 
as a result of the open string states ending on them. Conformal invariance implies that pinched surfaces, with infinitely long 
thin strips, have to be taken into account. In bosonic string theory, such surfaces can be resummed p^ . 



of the bosonic sigma-model, describing the interaction of string matter with D-particle defects. In section IIVI we 
discuss re-summation of higher world-sheet topologies and their implication on inducing stochastic fluctuations in 
the recoil velocity of the D-particle. The latter are expressed through a Langevin-type stochastic equation, which 
replaces the tree-level world-sheet renormalization group equation for the recoil velocity, interpreted as a sigma-model 
coupling. In section |Vl we solve this Langevin equation but also discuss "fuzzy" aspects of the string coupling in a full 
quantum-gravity setting of this string model of foam. For D-particle foam we thus arrive, in the sense of Tsallis, at a 
superstatistics description. We discuss the application of our results to a precise estimation of decoherence-induced 
CPT- violating effects in entangled states of mesons (the w-effect) in section IVI) where we estimate the effect in this 
model of (D-particle) space-time and compare it to previous estimates. Our conclusions and outlook are presented in 
section IVlIl Some technical aspects are reviewed in two Appendices. 



II. GENERIC FEATURES OF D-PARTICLE FOAM MODEL 



In this section we shall review briefly the basic features of the D-particle foam model, discussed in We will use 
some established results and constructs from string/brane theory 13, l5( , which we shall discuss briefly for the benefit of 
the non-expert reader. In particular, zero dimensional D-branes [6| occur (in bosonic and some supersymmetric string 
theories) and are also known as D-particles. Interactions in string theory are, as yet, not treated as systematically 
as in ordinary quantum field theory where a second quantised formalism is defined. The latter leads to the standard 
formulations by Schwinger and Feynman of perturbation series. When we consider stringy matter interacting with 
other matter or D-particles, the world lines traced out by point particles are replaced by two-dimensional world sheets. 
World sheets are the parameter space of the first quantised operators ( fermionic or bosonic) representing strings. In 
this way the first quantised string is represented by actually a two dimensional (world-sheet) quantum field theory. 
An important consistency requirement of this first quantised string theory is conformal invariance which determines 
the space-time dimension and/or structure. This symmetrypermits the representation of interactions through the 
construction of measures on inequivalent Riemann surfaces In and out states of stringy matter are represented 
by vertex insertions at the boundaries. The D-particles as solitonic states |5| in string theory do fiuctuate themselves 
quantum mechanically; this is described by stringy excitations, corresponding to open strings with their ends attached 
to the D-particles and higher dimensional D branes. In a first quantised (world-sheet) language, such fiuctuations are 
also described by Riemann surfaces of higher topology with appropriate Dirichlet boundary conditions (c.f. fig. [T]). 
The plethora of Feynman diagrams in higher order quantum field theory is replaced by a small set of world sheet 
diagrams classified by moduli which need to be summed or integrated over [l7|. The model of space-time foam, 
used in understanding the omega effect, is based on D-particles populating a bulk geometry between parallel D-brane 
worlds. The model is termed D-foam [4| (c.f. figure [2]), and our world is modelled as a three-brane moving in the bulk 
geometry; as a result, D-particles cross the brane world and appear for an observer on the brane as foamy structures 
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FIG. 2: Schematic representation of a D-foam. The figure indicates also the capture/recoil process of a string state by a 
D-particle defect for closed (left) and open (right) string states, in the presence of D-brane world. The presence of a D-brane 
is essential due to gauge flux conservation, since an isolated D-particle cannot exist. The intermediate composite state at 
t — 0, which has a life time within the stringy uncertainty time interval St, of the order of the string length, and is described 
by world-sheet logarithmic conformal field theory, is responsible for the distortion of the surrounding space time during the 
scattering, and subsequently leads to induced metrics depending on both coordinates and momenta of the string state. This 
results on modified dispersion relations for the open string propagation in such a situation Q, leading to non-trivial "optics" 
for this space time. 



Even at low energies E, such a foam may have observable consequences e.g. decoherence effects which may be of 

magnitude O ( j with n = 1,2 where Mp is the Planck mass or change in the usual Lorentz invariant dispersion 

relations. The study of D-brane dynamics has been made possible by Polchinski's realisation 5] that such solitonic 
string backgrounds can be described in a conformally invariant way in terms of world sheets with boundaries Q . On 
these boundaries Dirichlet boundary conditions for the collective target-space coordinates of the soliton are imposed 
[l9| . When low energy matter given by a closed string propagating in a. {d + l)-dimensional space-time collides with 
a very massive D-particle (0-brane) embedded in this space-time, the D-particle recoils as a result [l^ in a non- 
relativistic manner. We shall consider the simple case of bosonic stringy matter coupling to D-particles. Hence we 
can only discuss matters of principle and ignore issues of stability due to tachyons. However we should note that an 
open string model needs to incorporate for completeness, higher dimensional D-branes such as the D3 brane. This is 
due to the vectorial charge carried by the string owing to the Kalb-Ramond field. Higher dimensional D-branes (unlike 
D-particles) can carry the charge from the endpoints of open strings that are attached to them. For a closed bosonic 
string model the inclusion of such D-branes is not imperative (see figure [5]) although D-particle fluctuations would 
generally require them. The details of the higher dimensional branes are not essential for our analysis however. The 
current state of phenomenolgical modelling of the interactions of D-particle foam with stringy matter will be briefly 
summarised now. Since there are no rigid bodies in general relativity the recoil fluctuations of the brane and the 
effective stochastic back- reaction on space-time cannot be neglected. As we will discuss, D-particle recoil in the "tree 
approximation" i.e. in lowest order in the string coupling gs, is required to cancel infrared singularities in a higher 
order disc or Riemann sphere amplitude in open or closed string theory respectively; the recoil induces a non-trivial 
space-time metric. For e a positive infinitesimal, 0^ (t), the regularised step function is introduced, in terms of a 
contour integral 

ee(t) = ^ r ^e'-. (2.1) 

For closed strings colliding with a heavy (i.e. non-relativistic) D-particle, the metric has the form [2l| 

9ij = I 900 = -1, goi £ {eyi + Uit) Oe{t) , i = 1, . . . , d. (2.2) 

where the suffix denotes temporal (Liouville) components, Ui — (fci — ^2)^ is small, fci (fc2) is the momentum of the 
propagating closed-string state before (after) the recoil, are the spatial collective coordinates of the D- particle and 
is identified with the target Minkowski time t for t ^ after the collision. The latter requirement is consistent 
with e being infinitesimal. For our purposes the Liouville and Minkowski times can be identified. Now for large t, to 
leading order, 

_ _ u, Ap, 

got ~ u, = — oc g^— (2.3) 
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where Api is the momentum transfer during a colUsion and Ms is the string mass scale, < 1 is the string couphng, 
assumed weak, and the combination M^/qs is the D-particle mass, playing the role of the Quantum Gravity scale 
in this problem, i.e. the Planck mass; this formalism was used to establish a phenomenological model where the 
couplings Ui were taken to be stochastic and modeled by a gaussian process. The latter assumption was not based 
on analysis of the underlying string theory. The gaussian process represents a large universality class for stochastic 
processes and in this sense is an understandable assumption. Our purpose is to determine the justification of using 
gaussian (or other) distributions for the velocity recoil in the context of D-particle foam. 



III. MODULI AND VERTEX OPERATORS FOR D-PARTICLES: A COMPREHENSIVE REVIEW 



In order to understand the dynamics of D-particle foam it is imperative to consider D-particle quantisation and in 
particular D-particle recoil as a result of scattering off stringy matter. The issue of recoil is not fully understood and 
it is not our purpose here to delve into the subtleties of recoil and operators describing recoil. We will rather proceed 
on the basis of a proposal which has had some success in the past Q . The energy of a D-particle is independent 
of its position. Consequently in bosonic string theory there are 25 zero modes, the Dirichlet directions. Zero modes 
lead to infrared divergences in loops in a field theory setting where collective co-ordinates are used to isolate these 
infrared divergences. This is essentially due to the naivety of the perturbation series that is used and can be addressed 
using a coherent state formalism. The situation is similar but, in some ways, worse for string theory since the second 
quantised formalism for strings is more rudimentary. In string theory the use of first quantisation requires a sum 
over Riemann sheets with different moduli parameters. The formal transition from one surface to another of lower 
genus (within string perturbation theory ) has singularities associated with infrared divergences in the integration 
over moduli parameters. In the case of a disc D an incipient annulus E can be found by making two punctures and 
attaching a long thin strip (somewhat like the strap in a handbag and dubbed wormholes) . In the limit of vanishing 
width this wormhole, is represented by a region in moduli space, integration over which leads to an infrared divergence. 
A similar argument would apply to a Riemann sphere where two punctures would be connected by a long thin tube. 
For the disc the divergence will be due to propagation of zero mode open string state while for the Riemann sphere 
it would be an analogous closed string state along the wormhole. Let us examine this divergence explicitly. Consider 
the correlation function (V1V2V3 . . . Vn)^: for vertex operators [l| where S is the Riemann surface of an 

annulus; it can be deformed into a disc with a wormhole attached [l9|,[l8|. The set of vertex operators include 
necessarily any higher dimensional D-branes necessary to conserve string charge. However such vertex insertions 
clearly do not affect the infrared divergence caused by wormwholes. The correlation function can be expressed as 

{V1V2V3 ■■■Vn)^^J2 [ / / —q''--^ {^a (Sl) 0a (^2) Vi^s^S • • • K)z3 (3-1) 

^ J J J Q 

where si and S2 are the positions of the punctures on dD.The {0a} are a complete set of eigenstates of the Virasoro 
operator Lq with conformal weights hn |22l| and g is a Teichmuller parameter associated with the added thin strip. 
Clearly there is a potential divergence associated with its disappearance g — > and this corresponds to a long thin 
strip attached to the disc. For a static D-particle the string co-ordinates X d = {^^i ■ ■ • j X^^^ have Dirichlet 

boundary conditions while X'^ has Neumann boundary conditions (in the static gauge) 



oa 



^0 = ^X^{r,a) 



(3.2) 



where (r, cr) is a co-ordinisation of the worldsheet.The associated translational zero mode is given by 

0^(X,^) = ^a„Xl,e^'^^", (3.3) 

where dn denotes a derivative in the X^ Dirichlet direction, and is an element of the set {0a }• The conformal weight 
hi = 1 + a'uj"^ . The relevant part of the integral in (|3.ip is 

1 roo /.I 



dq dujq-'+'^- ^ dq^ — (3.4) 

J-oo Jo q (— \ogq) ' 
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and is divergent because of the behaviour of the integrand near g = 0. This can be regularized by putting a lower 
cut-off g > 5 — *■ in the integral. The correlation function on S can be computed to be (li^] 

{VlV2V:i . . . K)s = \0gS(^d„XD (Sl) .dnXo (S2) V1V2V3 . . . Vn)^ (3.5) 

where T is a cut-off for large (target) time. Division by it removes divergencies due to the integration over the 
world-sheet zero modes of the target time. These should not be confused with divergencies associated with pinched 
world-sheet surfaces, proportional to logS, that we are interested in here. These latter divergences cause conventional 
conformal invariance to fail. 

However, in one approach, it was argued sometime ago, that these divergences can be canceled if D-particles are 
allowed to recoil [isl . [20j . as a result of momentum conservation during their scattering with string states. In fact, 
in our D-particle foam model this is not a simple scattering process, as it involves capture and re-emission of 
the string state by the D-particle defect. In simple terms, this process involves splitting of strings by the defect. 
In a world-sheet (first quantization) framework, such processes are described by appropriate vertex operators, whose 
operator product expansion close on a (local) logarithmic algebra [H, |23| ■ The translational zero modes are associated 
with infinitesimal translations in the Dirichlet directions. Given that D-particles do not have any internal or rotational 
degrees of freedom, these modes should give us valuable information concerning recoil. Moreover (for w = 0) there is 
a degeneracy in the conformal weights between (j)^, d^(f)^ and the identity operator and also the conformal blocks in 
the corresponding algebra have logarithmic terms. 

To understand the formal structure of the world-sheet deformation operators pertinent to the recoil/capture process, 
we first notice that the world-sheet boundary operator Vd describing the excitations of a moving heavy DO-brane is 
given in the tree approximation by: 

/ {y^dnX' + u,X''dnX') = f Y, (XO) d^X' (3.6) 

JdD JdD 

where Ui and yi are the velocity and position of the D-particle respectively and (^°) = Vi + UiX'^. To describe 
the capture/recoil we need an operator which has non-zero matrix elements between different states of the D-particle 
and is turned on "abruptly" in target time. One way of doing this is to put a 6 {X'^^ , the Heavysidc function, in 
front of Vd which models an impulse whereby the D-particle starts moving at X° = 0. Using Gauss's theorem this 
impulsive Vd, denoted by V^™^, can be represented as 

v'^p = V / d^zdai [u,x°] e (X") = J][dT u,x°e (x^) dnX\ (3.7) 

»=1 -^D JdD 

Since X° is an operator it will be necessary to define 8 {X'^) as an operator using the contour integral 

(X°) ^-^ r ^^^e-^° with e ^ + . (3.8) 

Hence we can consider 

i?,(X") = Z?(XO; e) = X°Q, (^O) ^~ T ^ e'"^° . (3.9) 

J-oo (uj — ie) 

The introduction of the feature of impulse in the operator breaks conventional conformal symmetry, but a modified 
logarithmic conformal algebra holds. A generic logarithmic algebra in terms of operators C and V and the stress 
tensor T (z) (in complex tensor notation ) satisfies the operator product expansion 

T{z)C{w,w) ^^C{w,w) + ^P^^ + --- (3.10) 

(z ~ w) — w) 



T (z) V {w,W) ^ ^ V iw,w) + ^ C (w) + + • • • (3.11) 

(z — w) (z — w) [z — w) 



and 



(C(z,z)C(0,0)) -0 (3.12) 
1^1 

{V {z,z) V (0, 0)) ^ (log |z| + I)) (3.14) 

\z\ 
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where D is a constant. Since the conformal dimension of e^'^-^^' is ^ we find that 

T H (z) ^ - i^, (z) + -^-2 e©e M + • • • (3.15) 

2 (w — z) (w — z) 

and so a logarithmic conformal algebra structure arises if we define 

C,(X°) = C(X";e) = eee(X°) , (3.16) 

suppressing, for simplicity, the non-holomorphic piece. The above logarithmic conformal field theory structure is 
found with this identification. Similarly we find 

T{w)CAz)-^—^^CAz) + --- 
2[w — z) 

Consequently A for (z) and (z) is — A calculation (in a euclidean metric) for a disc of size L with a 
short-distance worldsheet cut-off a reveals that as e — > 

{C,{z)C,iO))^0(e^) (3.17) 



TT 



2 



(C,(z,z)i?,(0))~77A/T^ l-2£'log - (3.18) 



n TT /I 



2 



a 



D,{z,z)D,iO))^-J—[-~2\og- (3.19) 



where a = log | | . We consider e ^ 0-|- such that 



^2^^ 1 =0(1) , (3.20) 

where rj is the time signature and the right-hand side is kept fixed as the cutoff runs; it is then straightforward to see 
that ([XT7)) . dSIIHl), and (pi^ are consistent with (PTT^ . ([XT^ . and (PH)) . It is only under the condition ([X^ that 
the recoil operators and obey a closed logarithmic conformal algebra [20] : 

<C,iz)CM>-0 

< Ce{z)DAO) > - 1 

< De{z)De{0) > - -2rylog (3.21) 

The reader should notice that the full recoil operators, involving holomorphic pieces with the conformal- 

dimension-one entering (|3.7p ). obey the full logarithmic algebra p.l2p . (|3.13p . (|3.14p with conformal dimensions 

2 

A = 1 — From now on we shall adopt the euclidean signature ij = 1. 

We next remark that, at tree level in the string perturbation sense, the stringy sigma model (inclusive of the D- 
particle boundary term and other vertex operators) is a two dimensional renormalizable quantum field theory; hence 
for generic couplings it is possible to see how the couplings run in the renormalization group sense with changes 
in the short distance cut-off through the beta functions In the world-sheet renormalization group based on 
expansions in powers of the couplings, /3' has the form ( with no summation over the repeated indices) 

= y^9' + ■■■ (3.22) 

where yt is the anomalous dimension, which is related to the conformal dimension A^ by yi = Ai — S, with 5 the 
engineering dimension (for the holomorphic parts of vertex operators for the open string S = I). The ... in p.22p 
denote higher orders in g*. Consequently, in our case, we note that the (renormalised) D-particle recoil velocities 
constitute such a-model couplings, and to lowest order in the renormalised coupling Ui the corresponding /3 function 
satisfies 

-u\ (3.23) 



dlogA 2 

where A is a (covariant) world-sheet renormalization-group scale. In our notation, we identify the logarithm of this 
scale with a = log | | , satisfying p.20p . 
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An important comment is now in order concerning the interpretation of the flow of this world-sheet renormalization 
group scale as a target-time flow. The target time t is identified through t — 2 log A. For completeness we recapitulate 
the arguments of [20| leading to such a conclusion. Let one make a scale transformation on the size of the world-sheet 

L^L' = e'/^L (3.24) 

which is a finite-size scaling (the only one which has physical sense for the open string world-sheet). Because of the 
relation between e and L (|3.20p this transformation will induce a change in e 

(note that if e is infinitesimally small, so is e' for any finite t). From the scale dependence of the correlation functions 
p.21[) that Ce and transform as: 

C, ^ C,, = C, (3.26) 

From this transformation one can then see that the coupling constants in front of and in the recoil operator 
(13. 6p . i.e. the velocities Ui and spatial collective coordinates Ui of the brane, must transform like: 

Ui , yi —>yi+ Uit (3.27) 

This transformation is nothing other but the Galilean transformation for the heavy D-particles and thus it demon- 
strates that the finite size scaling parameter i, entering ()3.24p . plays the role of target time, on account of (j3.20p . 
Notice that (|3.27p is derived upon using (|3.2ip . that is in the limit where e ^ 0. This will become important later 
on, where we shall discuss (stochastic) relaxation phenomena in our recoiling D-particle. 

Thus, in the presence of recoil a world-sheet scale transformation leads to an evolution of the Z?-brane in target 
space, and from now on we identify the world-sheet renormalization group scale with the target time t. In this sense, 
equation (j3.23p is an evolution equation in target time. 

However, this equation does not capture quantum-fluctuation aspects of about its classical trajectory with time 
Ui{t). Going to higher orders in perturbation theory of the quantum field theory at fixed genus does not qualitatively 
alter the situation in the sense that the equation remains deterministic. In the next section we shall consider the 
effect of string perturbation theory where higher genus surfaces are considered and re-summed in some appropriate 
limits that we shall discuss in detail. 



IV. STRING PERTURBATION THEORY AND IMPLICATION FOR RECOIL VELOCITY 



It is not possible to exactly sum up higher orders in string perturbation theory. We have seen that infrared 
singularities in the integration over the moduli of the Riemann surface (representing the world sheet) in the wormhole 
limit are related to the recoil operators for the D-particle. The wormhole construction [24] is a way of constructing 
higher genus surfaces from lower genus ones. Since it will be relevant to us later, we should note that gs the string 
coupling is given by 

c;,=eW (4.1) 

where $ is the spin zero dilaton mode which is part of the massless string multiplet. Here (. . .) denotes the string 
path integral / DX e^' $ where is the string ct- model action in the presence of string backgrounds such as the 
dilaton and the Kalb-Ramond modes. In particular the a model deformation due to the dilaton has the form 

— / dadr y^<i>{X)R^^'> {t,(j) (4.2) 

on a worldsheet Riemann surface S where jafs is the induced metric on the worldsheet, 7 — \detjai3\ and i?^^' is the 
associated Ricci curvature scalar. Now the Euler characterisitic x of S is given by 

X = ^ / dadT^R^^^ = 2 (1 - ff) (4.3) 
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where g is the genus and is an integer valued invariant. If we spUt the dilaton into a classical (worldsheet co-ordinate 
independent) part ($) and a quantum part = : $ : , where : . . . : denotes appropriate normal ordering, we can 
write <f> = ($) + The cr-model partition function Z can be written as a sum over genera 



Z = Y, / d7„;3dXe-^-*-x<*>-i^-'^'^''^V7^fl<^'(-.'^) (4.4) 
X 

= J29s^ f /rf7a/3d^e-^-'-3^/B'^'^''^^^«'"(^^'^) (4.5) 



X 

where Srest denotes a u-model action involving the rest of the background deformations except the dilaton. For the 
moment we will assume that the theory is such that a potential is generated for $ which suppresses the fluctuations 
represented by (p. In general we would have to consider — e* which would then make the string coupling a field. 

The summation over genera cannot be performed exactly. We will follow an approach using a mechanism due 
to Fischler and Susskind based on a dilute gas of wormholes (proposed originally by Coleman within the 

context of euclidean quantum gravity [131 )■ This results in the structure of recoil (in lowest order) being modified by 
generating a gaussian distribution for the recoil velocity u\ 

We present a detailed review of the pertinent formalism in Appendix A. For our purposes in this section we note 
that, in the case of mixed logarithmic states, the pinched topologies are characterized by divergences of a double 
logarithmic type (c.f. (|7.7p in Appendix A) which arise from the form of the string propagator (c.f. (|7.ip in Appendix 

A) in the presence of generic logarithmic operators C and D, J dq q"^'^^ {C,D\ ^'-'S'?^ \C,D) . 

As shown in [13], the mixing between C and D states along degenerate handles leads formally to di- 
vergent string propagators in physical amplitudes, whose integrations have leading divergences of the form 
logg/ d?z D{z\e) f cPz' C{z'^ f (Pz D{z;e) J cfz' C{z';e) .As explained in [3, and reviewed in Appendix 
A of the current manuscript, these (log (5)^ divergences can be cancelled by imposing momentum conservation in the 
scattering process of the light string states off the D-particle background. 

We note at this stage that, as mentioned earlier, isolated D-particles do not exist, as a result of their gauge flux 
conservation requirement. The physically correct way to formulate, therefore, the problem, is to consider groups 
of TV, say, D-particles, which interact among themselves with flux-carrying stretched strings. The above analysis 
remains intact (in the sense of generalising straightforwardly) when more than one D-particle is present in a region 
with typical dimensions smaller than the string length sometimes known as the fat brane [isj . There are, of course, 
technical differences in the sense that a non-abelian structure arises and the Y have matrix labels. All qualitative 
features, however, are preserved concerning the quantisation of the D-particle background moduli. For the remainder 
of this section we shall, therefore, formulate our arguments within this rigorous multi-D-particle picture. 

The cancelation of leading divergences of the genus expansion in the non-abelian case of a group of N D-particles, 
is demonstrated explicitly in appendix A. It is shown there that this renormalization requires that the change in 
(renormalized) velocity of the due to the recoil from the scattering of string states be 



1 / \ rIY"''' 

- (k, + k2) ^ , a,b = l...N (4.6) 

L \ J i at 



Ijab ^ ■ ■ 
M. 



where fci.2 are the initial and final momenta in the scattering process and Mo — 1/^/a' g^ is the BPS mass of the 
string soUton [1, [H], and gs < 1 is the physical (weak) string coupling. In (|4.6p . the fci,2 are true physical momenta 
so that Mu represents the actual BPS mass of the D-particles. This means that, to leading order, the constituent 
D-particles in a group of N of them, say, move parallel to one another with a common velocity and there are no 
interactions among them. Thus the leading recoil effects imply a commutative structure and the "fat brane" of the 
group of D-particles behaves as a single D-particle (with a single average collective coordinate of its center of mass). 
In such a limit one may replace J/f'' by Ui (c.f. previous section), describing the collective recoil velocity of the fat 
brane. This should be understood throughout this work. 

In addition to this divergence, there are sub-leading log (5 singularities, corresponding to the diagonal terms 
/ (Pz D{z;e) J (Pz' D{z';e) and J(Pz C{z;e) J cPz' C{z';e). These latter terms are the ones we should concen- 
trate upon for the purposes of deriving the quantum fluctuations of the collective D-particle coordinates. It is these 
sub-leading divergences in the genus expansion which lead to interactions between the constituent D-branes and pro- 
vide the appropriate noncommutative quantum extension of the leading dynamics (|4.6p . The reader should recall that 
these (sub-leading) divergences also showed up in the much simpler case of perpetual Galilean motion of D-branes 
discussed in [l^ (c.f. p.Sp ). as a result of the translational symmetries zero mode contributions. 

In the weak-coupling case, we can truncate the genus expansion to a sum over pinched annuli (fig. [3] in Appendix 
A) . This truncation corresponds to a semi-classical approximation to the full quantum string theory in which we treat 
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the D-particles as heavy non-relativistic objects in target space. Then the dominant contributions to the sum are 
given by the log 5 modular divergences described above, and the effects of the dilute gas of wormholes on the disc are 
to exponentiate the bilocal operator (|7.ip of Appendix A, describing string propagation in a pinched annulus. Thus, 
in the pinched approximation, the genus expansion of the bosonic cr-model leads to an effective change in the matrix 
(7-model action by [18|] 



AS*: 



9. 



2 i°g^ 2. 



duj duj' 



as JdT. 



(4.7) 



where w, w' are Fourier variables, defined appropriately in Appendix A, and , i, j = C, Z) is a metric in the theory 
space of strings, introduced by Zamolodchikov [25|. 

The bilocal actio n ( 1 4. 7p can be cast into the form of a local worldsheet effective action by using standard tricks of 
wormhole calculus [24l| and rewriting it as a functional Gaussian integral [l8l | 



exp 



gs E / /^f (c^) f V:,{a 



(4.8) 



where pf^{u!) are stochastic coupling constants of the worldsheet matrix ct- model, which express quantum fluctuations 
of the corresponding background fields in target space, as a consequence of genus re-summation. Thus the effect of 
the resummation over pinched genera is to induce quantum fluctuations of the collective D-brane background, leading 
to a set of effective quantum coordinates 



viewed as position operators in a co-moving target space frame. 

Thus we find that the genus expansion in the pinched approximation is [Tsj 



E Zf' [A] 0, If [dp] p[p] W [dH- A-4^p\) 



(4.9) 



(4.10) 



where the sum is over all pinched genera of infinitesimal pinching size, and 



p[p] cx exp 



1 

2r2 



E ds ds' pf {X'{s)) G%^,{s, s') pf (X0(.')) 

u . ^ Jo 



a,b.c,d 



is a (appropriately normalized) functional Gaussian distribution on moduli space of width 



r = gs \/log6 



(4.11) 



(4.12) 



In (j4.10p we have normalized the functional Haar integration measure [dp] appropriately. 

We see therefore that the diagonal sub-leading logarithmic divergences in the modular cutoff scale S, associated with 
degenerate strips in the genus expansion of the matrix a-model, can be treated by absorbing these scaling violations 
into the width F of the probablity distribution characterizing the quantum fluctuations of the (classical) D-brane 
configurations Y"''' {X'^ {s)) . In this way the interpolation among families of D-brane field theories corresponds to a 
quantization of the worldsheet renormalization group flows. Note that the worldsheet wormhole parameters, being 
functions on the moduli space of recoil deformations, can be decomposed as 



pfiX^s)) = hm {[pc]fC{X°;e) + i?(X"; e)) 



(4.13) 



The fields pc,D are then renormalized in the same way as the D-brane couplings, so that the corresponding renormal- 
ized wormhole parameters generate the same type of (Galilean) /3- function equations (|3.23[) . 

According to the standard Fischler-Susskind mechanism for canceling string loop divergences i 28:] , modular infinities 
should be identified with worldsheet divergences at lower genera. Thus the strip divergence log S should be associated 
with a worldsheet ultraviolet cutoff scale log A, which in turn is identified with the target time as described earlier. 
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We may in effect take 5 independent from A, in which case we can first let e — > 0+ in the above and then take the 
hmit S 0. Interpreting log (5 in this way as a renormalization group time parameter (interpolating among D-brane 
field theories), the time dependence of the renormalized width (|4.12p expresses the usual properties of the distribution 
function describing the time evolution of a wavepacket in moduli space. The inducing of a statistical Gaussian spread 
of the D-branc couplings is the essence of the quantization procedure. 

A final remark is in order. From the form ()3.9p and (j3.16|) of the recoil operators, it is evident that the dominant 
contributions in the limit e — > 0+, we consider here, come from the Z?-deformations, pertaining to the recoil velocity 
of the D-particle (or, better, the center of mass velocity of a group of D-particles, as discussed above). From now 
on, therefore, we restrict our attention to the distribution functions of such recoil velocities: 

1 



p(«)^_e-^^, r = g,yi^, (4.14) 

where u denotes the classical recoil velocity. Notice that, upon invoking [l^ the Fischler-Susskind mechanism [2^ 
for the absorption of the modular infinities to lower-genus (disc) world-sheet surfaces, we may identify logS with the 
target time: 

\ogS = t , (4.15) 

where this identification should be understood as being implemented at the end of the computation. To be precise, 
as explained in [l^, the correct form of (|4.15p would be: \og5 = g} t, with x > an exponent that can only be 
determined phenomenologically in the approach of |l3|, by comparing the space-time uncertainty principles, derived 
in this approach of re-summing world-sheet genera, with the ones within standard string/brane theory. In fact, in 
our approach of re-summing world-sheet pinched surfaces [Tsj . one obtains for the spatial and temporal variances: 
AF°°At > ^Ja/c/, which implies that the standard string-theory result [13], independent of the string coupHng, is 
obtained for x = 0. This is the case we shall consider here, which leads to the identification (|4.15p . However, in the 
modern approach of D-brane theories, one can adjust the uncertainty relations in order to probe minimal distances 
below the string length, which is achieved by the choice [3, e.g. x = 2/3, reproducing the characteristic minimal 
length probed by D-particles [Slj . In our case, where, as we shall discuss in the next subsection, the coupling constant 
of the string may itself fluctuate, it is the mean value of gs that enters in such relations. This issue is not relevant if 
we stay within the x = case, which we do in this article. 

We next remark that the nature of the Gaussian correlation is assumed to be delta correlated in time. The Langevin 
equation [IHl implied by ()4.9p replaces (|3.23|) and can be written as 

^^-lu'-^^W^at) (4.16) 
at U ^rin/i 



where t = and ^ (i) represents white noise. This equation is valid for large t. From the above analysis it is known 
that [18] (c.f. Appendix A)) that to O (g^) the correlation for ^ (t) is independent, and for time scales of interest, 
is correlated like white noise ; hence the correlation of ^ (i) has the form: 

(<) e (*')> = '5 (i-o- (4.17) 

Since the vectorial nature of is not crucial for our analysis we will suppress it and consider the single variable u. 

We should stress that this equation is valid for large t which is required since e is small. Hence the apparent 
singularity in Ean (j4.16l ) at i = is not relevant and so we can empirically regularise this singularity by changing i 
to for some to > 0; to is the order of the capture time of the (j) meson by the D-particle The stochastic Langevin 
equation ()4.16p . describes relaxation aspects of the recoiling D-particle with equilibrium being reached only as e — )■ 
(or t oo). 

The reader should notice that in the limit the system reaches equilibrium with a constant in time velocity. It is 
only in this limit that the Galilean transformation (|3.27p applies, as already discussed there. We now proceed to a 
solution of this Langevin equation and a discussion on the pertinent physical consequences for a statistical population 
of quantum-fluctuating D-particles. 



V. SOLUTION OF LANGEVIN EQUATIONS AND FLUCTUATING STRING COUPLING 



In the recent modeling of the omega effect [9|, the recoil velocity of the D-particles u has been taken as a classical 
stochastic variable. The D-particle fluctuations which are described by (|4.16p will be superimposed on this stochas- 
ticity. Eqn. (|4.16p is particulary simple equation in the sense that the drift and diffusion terms are independent 
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of u. By making a change of variable it is easy to eliminate the drift term and the resulting equation can then be 
interpreted in terms of a Wiener process p6| . Let us consider the auxiliary equation 



1 



y, 



dr 4(t + tn) 
which just deals with the drift part of Eqn. (|4.16p . It has a solution 

y{t)=y{t^)T{t) 

where 



(5.1) 



T [t) = exp 



dt' 



1 



4 7o t' + to 



t+Jo 

to 



and to is a time much smaller than t. We now define U (t) = u {t) T {t) ^ and readily find that 



dU_ 
'dt 



/2a' 



■.t^/^T{t)-'^{t). 



(5.2) 



(5.3) 



This describes purely diffusive motion and is thus related to the Wiener process; equivalently we can consider the 
associated probability distribution p {U, t) which satisfies the Fokker-Planck equation 



d 



1 



to ' 



If at i = consider a D-particle velocity recoil uq so that 

p([/,0) = <5({7 

The Eqn. (|5.4p can be solved to give 



v{u,t) 



I 15a' 1 
27r?7 (t) 5, 



■ exp 



15a'{U -uof 



where 



2tl 



i{t + tof 



1 



to 



5t^ {t + to)^ 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



If the D-particle is typically interacting with matter on time scales of to, then the effect of a large number of such 
collisions can be calculated by performing an ensemble average over a distribution of uo- A distribution for uo that 
has been used in modelling is a gaussian with zero mean and variance a. This is readily seen to lead to an averaged 
distribution D-particle velocity recoil distribution <^ p{u\gs) ^ where 



< p (u ) > 



(5.8) 



15a' 



27r {glri {t) + Iba'a^ 



■ exp 



15a' 



2 iSsV (t) + I5a'a^ 



We have used a notation for p which emphasises that it is conditional on gs having a fixed value. The interaction time 
includes both the time for capture and re-emission of the string by the D-particle, as well as the time interval until 
the next capture, during string propagation. In a generic situation, this time could be much larger than the capture 
time, especially in dilute gases of D-particles, which include less than one D-particle per string (a^^^) volume. Indeed, 
as discussed in detail in |32| , using generic properties of strings consistent with the space-time uncertainties [30] , the 
capture and re-emission time to, involves the growth of a stretched string between the string state and the D-brane 
world (c.f. fig.[2|) and is found proportional to the incident string energy p'^: 



to 



a'p° < x/c/ 



(5.9) 
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We shall use this result in section IVIl where we estimate the strength of the w-effect in the initial entangled state of 
two mesons (jl.2p . after the 0- meson decay in the presence of D-particles. In such a situation, the interaction time is 
essentially the capture time to- 

We will now examine how the above results are modified when gs fluctuates. Such issues cannot currently be 
treated with any level of rigour since they embody issues of string vacua and backgrounds. Hence we will take a 
somewhat phenomenological stance and consider the effect of a class of stochastic fluctuations for g^^ (= c^^^*^) 
which are varying on a timescale which is slow compared to the drift time-scale. This is the arena of superstatisitcs 
[l^ . Following the analysis of Beck [27j] (c.f. Appendix B) we can make the following fairly general classical ansatz 



for 



compatible with its positivity viz. 



9t 



1=1 



(5.10) 



where the Xi are n independent gaussian variables of zero mean and variance CTq . The probability distribution p for 



is X with n degrees of freedom, i.e. 



".9o 



i/2 



X I exp 



"5o 



(5.11) 



Here 



-ir and the variance var 

9o 



If go is held constant the variance can be made small for large 



2 

n. This is a particular form of the Gamma distribution. However other choices for p, such as lognormal and F 
distributions have also been considered [ij] in contexts such as turbulence. We note that 



e-""- ''/3^-iexp 
We now calculate the probability p is for u as 

p{u) = 



n(3 
Wo 



d(3 = 



n 

Wo 



< p (u ) > . 



(5.12) 



(5.13) 



The choice of p in (|5.1ip leads to canonical form of distribution for non-extensive statistics jl5|,[27[ . Other choices 
mentioned above give different superstatistics. These wider classes of statistics present opportunities in interpreting 
the data in neutrino physics (see e.g. ^^)- For large n it is straightforward to show that 



15q;' 



T:{glil{t) + lba'n^) 



p{u) 



On writing q = 1 + and g (t) = {3 — q) la 



i/2 



(n\ ( ISa'gp 
y2) \i-,(t)gl + lZa' 



ngl 



(5.14) 



15ck' 



p{u) 



155 (i) 



we find the canonical form for Tsallis statistics [15 1 

(g-i)^i' 



(5.15) 



with q the non-extensivity parameter. Hence the stochasticity in the recoil velocity, when stringy matter is captured 
by the D-particle, leads to a deviation from nonextensive Tsallis statistics. The deviation is suppressed as the 
interaction time increases. For large n we have weak non-extensivity and the fluctuations for ^ are small. 

The quantity of interest which is an important input for our estimate of the omega effect is the variance of p (u) : 



var(u) 



T] it) gl + Iba'a'^ 



15a' 1 



(5.16) 



For reasons mentioned above, we can make the plausible assumption that the initial state of the neutral K meson 
pairs is governed by the variance a.t t ^ to ^ a'p^. In such a case, from (15. 7p we obtain that rj{t ^ to) — 2(1 + \/2)to ^ 
2{1 + V2){a'p°f, and hence from ([5T6l) the variance over such time scales becomes of order: 



Uar(u)non-cxt {t to) 



2{l + V2)[gUV ^p' 
15(1 



.0\2 



goV(pY(l + - + ---)+0(^') 



(5.17) 
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where on the right hand side of the above equation we only gave an order of magnitude estimate, assuming that n is 
large. The terms of order have not been written explicitly. Indeed, in most models of quantum gravity, a natural 
assumption would be that < g^a' (p'^Y , which is a natural assumption to make for a dispersion due to (quantum) 
fluctuations of the recoil velocity of heavy D-particles of average mass Ms/go = ^/igoVa'), where {^/g1) — l/.9o- As 
we shall discuss in the next section, therefore, such dispersion terms do not lead to dominant contributions to the 
w-effect estimates which are of primary interest to us here. These considerations lead to natural estimates for the 
parameter C, in the magnitude of the omega effect that we have discussed in our earlier work [ol] , and we now proceed 
to examine. 



VI. DECOHERENCE AND ENTANGLED STATES: (.j-EFFECT REVISITED 

We shall use a low-energy quantum-mechanical approach for the dynamics of the neutral mesons, which is sufficient 
for a discussion of the w-effect in non-relativistic systems of entangled mesons, such as Kaons in a (/)-factory Q. For 
a description of the stochastically fluctuating space-time effects, we shall make use of the above-derived string theory 
effects, in particular the recoil- velocity dispersion (|5.17p . including the effects of non-extensive statistics (jS.lSp . due 
to fluctuations of the string coupling gs ■ 

Following 0, where we refer the interested reader for details, we consider the following interaction Hamiltonian, 
which expresses the effective low-energy interaction of the meson states with the D-particle foam space-time back- 
ground in the model of [4;] : 

where k indicates the appropriate momentum operator, with eigenvalue k (along the direction of motion), when 
acting on momentum eigenstates, i.e. k |±fc, |) = ±fc |fc, ]) together with the corresponding relation for |. The arrows 
indicate the appropriate meson "flavours" |9]. The induced metric g^^, on the other hand, is such that: 

5"" = (-l + r4)l 

50i=g,io = rol + riCTi+r2CT2+r3cr3 (6.2) 
5" = (1 + ^5)1 

where 1 , is the identity and a.i ,i = 1, 2, 3 are the Pauli matrices. 

The target space metric state, which is close to being flat, can be represented schematically as a density matrix 

Pgrav= I d\ f{r,,)\g{r^)){g{r^)\. (6.3) 

The parameters {fj, — 0, ... ,5) are stochastic with a gaussian distribution / (r^ ) characterised by the averages 

(r^) = 0, (vr,) = A^5^, . (6.4) 

The fluctuations experienced by the two entangled neutral mesons will be assumed to be independent. 

The above parametrisation has been taken for simplicity and we will also consider motion to be in the x- direction 
which is natural since the meson pair moves coUinearly in the Center-of-Mass frame. 

Space-time deformations of the form (|6.2p and the associated Hamiltonians (|6.ip have been derived in the context 
of conformal field theory in earlier works by one of the authors and collaborators [1, [s^ and details on the relevant 
derivations will not be given here. We only mention that the variable ri in particular, expresses a momentum transfer 
during the interaction of the (string) matter state with the D-particle defect. In this sense, the off-diagonal metric 
component 501 can be represented as 

501 ~ Ul (6.5) 

where Ui = gs^ expresses the momentum transfer along the direction of motion of the matter string (taken here to 
be the x direction). In the above equation, Ms/ gs is the mass of the D-particle, which for weakly coupled strings with 
coupling gs is larger than the string mass scale Ms. In order to address oscillation phenomena, induced by D-particles, 
the fluctuations of each component of the metric tensor are taken in Q to have a 2 x 2 ("flavour") structure, as in (|6.2p . 
and not the simple structure (j6.5p considered in earlier works. In the case of neutral Kaons, which we concentrate on 
for concreteness in this section, we use the following notation for the "flavours" : 



\Kl) = IT) , \Ks) = U) 
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which represent the two physical eigenstates, with masses mi = m^, m2 = msj with 

Am^TriL-ms ^ 3.48 x 10"^^ GeV . (6.6) 

In this way, the stochastic variables f (|6.4p ) in (|6.2p . are linked with the fluctuations of the D-particle recoil velocity, 
by representing the latter as: 

(6.7) 

upon the above-mentioned technicality of considering flavour changes in addition to the momentum transfer. Specif- 
ically, \n\=0{gsrk/Ms). 

In this sense, the detailed discussion in the previous session on the stochastic fluctuations of the recoil velocity 
about a zero average value, translates into rewriting (|6.4p with variances (c.f. (|5.17p ) 



where we considered the capture time to a'p^ , with the energy of the probe, as spanning the essential interaction 
time with the D-particle of the initial entangled meson state. In this way we extrapolate the result (j5.8p to times 
smaller than Vc/. This is acceptable, as long as such times are finite. From a conformal fleld theory point of view, this 
means that we consider the world-sheet scaling parameter 1/e^ ~ \n{L/af' ^ to \/a' for probe energies <C Ms. 
We next note that the Hamiltonian interaction terms 

Hi = - (nai + r2a2)k (6.9) 

are the leading order contribution in the small parameters in the Hamiltonian H (|6.ip . since the corresponding 
variances y^A^ are small. The term (16. 9p . has been used in [9| as a perturbation in the framework of non-degenerate 
perturbation theory, in order to derive the "gravitationally-dressed" initial entangled meson states, immediately after 
the (j) decay. The result is: 



where 



+ \k, \-k, if^ - /3(2)) + |fc, T)(^^ \-k, ^f'> - ad); 

+ |fc,i)« |-fc,T)''' |fc,T)^'^ \-kAf^ (6.10) 



a 



E2 — El El 



P^^^ ^ ' ' , ^-1,2 (6.11) 



where the index (i) runs over meson species ("flavours") (1 Kl, 2 Ks). The reader should notice that the terms 
proportional to (a*^^) — a (1') and (/?(!) - /3(2)) in (l6l0l) generate w-like effects. We concentrate here for brevity and 
concreteness in the strangeness conserving case of the cj-effect in the initial decay of the 4> meson Q , which corresponds 
to ri oc di2. We should mention, however, that in general quantum gravity does not have to conserve this quantum 
number, and in fact strangeness- violating w-like terms are generated in this problem through time evolution Q. 

We next remark that on averaging the density matrix (c.f. (|6.3p ) over the random variables r.^, which are treated as 
independent variables between the two meson particles of the initial state (|6.10p , we observe that only terms of order 
Icjp will survive, with the order of |wp being 



^(l),(2) 



(6.12) 

for the physically interesting case of non-relativistic Kaons in (j) factories, in which the momenta are of order of the 
rest energies. The notation (2) (...) above indicates that one considers the sum of the variances A2 over the two 
meson states 1, 2 as defined above. 

The variances in our model of D-foam, which are due to quantum fluctuations of the recoil velocity variables about 
the zero average (dictated by the imposed requirement on Lorentz invariance of the string vacuum) are given by (j6.8p , 
with p*^ rrii the energy of the corresponding individual (non-relativistic) meson state (i), i = 1,2, in the initial 
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entangled state (|6.10p . It is important to notice that, on taking the sum of the variance A2 over the mesons (1) 
and (2), the terms proportional to the dispersion in the initial recoil velocity uq Gaussian distribution in (|5.17p 
give a contribution of order SOcr^, since is assumed universal among particle species. This is a parameter that 
depends on the details of the foam. As already mentioned in the previous section, one may assume models in which 
<C {goVa' p'^Y ■ In this sense, one is left with the contributions from the first term of the right-hand-side of (|5.17p . 
and thus we obtain the following estimate for the square of the amplitude of the (complex) w-parameter: 



1^1- go ' \,2 7 ^ 1 + ' Mslga = Mp, (6.13) 

(mi — m2) \ n J 

where Mp — Ms/go is the (average) D-particle mass, as already mentioned, representing the (average) quantum 
gravity scale, which may be taken to be the four-dimensional Planck scale. In the modern version of string theory, 
Mg is arbitrary and can be as low as a few TeV, but in order to have phenomenologically correct string models with 
large extra dimensions one also has to have in such cases very weak string couplings go, such that even in such cases 
of low Ms, the D-particle mass Ms/ go is always close to the Planck scale 10^^ GeV. But of course one has to keep an 
open mind about ways out of this pattern, especially in view of the string landscape. 

The result (|6.13p . implies, for neutral Kaons in a factory, for which (|6.6p is valid, the estimate \uj\ = O (lO^^), 
which, in the sensitive 77"' bi-pion decay channel, leads to effects enhanced by three orders of magnitude, as a result 
of the fact that the |w| effect always appears in the corresponding observables ,8] in the form |, and the 

CP-violating parameter |?7^ | ~ 10~^. At present, this value is still some two orders of magnitude away from current 
bounds of the w-effect by the KLOE collaboration at DA$NE giving \lo\ < 10~^ but it is within the projected 
sensitivity of the proposed upgrades. 

The above estimate is valid for non-relativistic meson states, where each meson has been treated as a structureless 
entity when considering its interaction with the D-foam. The inclusion of details of its strongly-interacting substructure 
may affect the results. To consider relativistic meson states, the non-relativistic quantum mechanics formalism leading 
to (|6.13p in the Kaon systems should strictly be replaced by an appropriate relativistic treatment. The major difference 

in this case is the form of the Hamiltonian, which stems from the expansion of the Dirac Hamiltonian for momenta 

2 

fc ^ mi, rrii the masses. The quantities Ei ^ k -\- « = 1, 2 (due to momentum conservation, assumed on average), 
and the capture times t,, ~ a'Ei. Nevertheless, if one extrapolates naively the above results for the |wp in such a case, 
one arrives at the conclusion that for the validity of leading-order perturbation theory for the interaction with the 
D-foam, one should consider momenta of order fc 30 GeV. Otherwise higher-order corrections become important. 
This completes our discussion on the estimates of the w-effect in the initial entangled state of two mesons in a meson 
factory. As discussed in [3], w-like terms can also be generated due to the time evolution. One can apply similar 
estimates for this case too. We shall not do so in this work. 



VII. CONCLUSIONS AND OUTLOOK 



This work examines the role of quantum string fluctuations (which can give rise to a non-commutative space-time 
geometry at string scales) on the velocity distribution of D-particles, within a specific kind of foam in string theory. 
Our Gaussian modeling of the recoil velocity is found to be robust to these fluctuations. 

In this way we have managed to give a rather rigorous estimation (modulo strong interaction effects) of the cj-effect 
in entangled states of mesons, which has been compared to previous naive estimates, based on dimensional analysis. 
The effect is smaller by roughly two orders of magnitude from the current upper bounds set by the KLOE collaboration 
at DA<I>NE. This makes the detectability of the effect in future meson factories a major experimental challenge. 

Admittedly, our approach in this paper is based on bosonic string theory which is not the most relevant phenomeno- 
logically. World-sheet supersymmetric strings do not lead to a closed-form resummation of the leading divergencies 
of the pinched surfaces, since the latter cancel out [21I, [s^, and the remaining terms are hard to cast in a closed 
form. However, despite this apparent technical difficulty, the general conclusions drawn from the current work, as 
far as fuzzyness of the string coupling and the target space time are concerned, are likely to be robust, since they 
depend on the form of vertex operators for (recoil) zero modes of D-particles (provided of course the theories are 
restricted to those admitting D-particles). In this sense, our approach here opens up the possibility of extending the 
analyses of D-particle foam to non-extensive statistics. The latter as a class is ubiquitous but, as far as we are aware, 
the possibility of non-extensive statistics in D-particle recoil has not hitherto been raised. Moreover, the mechanism 
makes contact with some non-trivial features of string theory. 

The role of such non-extensivity on matter propagation in D-particle foam, matter number distribution functions 
and (supercritical) string cosmology ^35^ will be addressed in future publications. This is an important aspect of the 
formalism discussed here, since it may have a profound influence on the dark matter (and dark energy) distributions 
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in a Universe with D-particle foam, which may have phenomenological consequences, as far as constraints on, say, 
supersymmetric particle physics models are concerned. Indeed, as we have discussed above, the interaction of D- 
particles with matter leads to local distortions of the neighboring space- time (|2.3p . which depend on both the string 
coupling gs, through the D-particle masses Ms/gs, and the recoil velocity of the D-particle (i.e. the momentum transfer 
of the string matter) Ui, which stochastically fluctuates upon summing up higher-genus world-sheet topologies. In 
view of our discussion in this work, both these quantities can be "fuzzy" , leading to stochastic fluctuations on the 
space-time metric, on which string matter lives. These fluctuations are up and above any statistical fluctuations in 
populations of D-particles that characterise the D-particle foam models. 

The presence of such fluctuations affect important cosmological quantities that are directly relevant to the (supercrit- 
ical) Universe budget, such as thermal supersymmetric dark matter relic densities, through appropriate modifications 
of the relevant thermodynamic equations. Hence, the relevant astro-particle physics constraints on supersymmetric 
models [sl] are also modified. As we have noted above, however, the issue as to whether the fuzzyness of the string 
coupling and consequently the non-extensivity of the D-particle foam can lead to observable signatures in Cosmology 
or astro-particle physics in general, remains to be seen. We hope to be able to report in a more detailed form on these 
phenomenological issues in the near future. 
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Appendix A: Recoil and Leading Divergences in the World-Sheet Genus Expansion 

In this appendix we shall show how the leading (log (5)^ modular divergences which appear in (|7.8p can be removed 
by invoking an appropriate Ward identity for the fundamental string fields of the matrix cr-model. As we shall show, 
this is equivalent to imposing momentum conservation for scattering processes in the matrix D-brane background. 
We shall be brief in our discussion, restricting oursleves on the main results, of relevance to our discussion in this 
work. The material in this section is taken from [isj , where we refer the interested reader for further details. 

As shown in that reference, conformal invariance requires absorbing such singularities into renormalized quantities 
at lower genera, leading to a generalized version of the Fischler-Susskind mechanism [28]. Such degenerate Riemann 
surfaces involve a string propagator over thin long worldsheet strips of thickness S ^ that are attached to a disc. 
These strips can be thought of as two-dimensional quantum gravity wormholes. Consider first the resummation of one- 
loop worldsheets, i.e. those with an annular topology, in the pinched approximation (c.f. fig. [3]). String propagation 
on such a worldsheet can be described formally by adding bilocal worldsheet operators B, which in the present case 
are defined by 

B{u;,u;')^Y. f f '\ , V^A^^^') (7-1) 

a,b.c.d J 9^ •'9^' ^0 ~ ^ 

where Gij denotes the Zamolodchikov metric in string theory space, that is the two-point correlation function of the 
recoil vertex operators: 

/ U„\(a;;^)^ Xi^^Y, f dsu e"^'^''" '^^^'^ Q{X\sk)-,e)Usk - e)Usk) ^X\su) . (7.2) 
and Lq denotes the usual Virasoro generator. The variable a; is a Fourier variable, appearing in the Fourier transform 

/•oo 

Y^^^{uj)^ Mm / dte'^'Yt\X°-e) (7.3) 

of the collective coordinates 1^"''(A°;£) of the (group of) D-particles (isl. [20j 

y,"''(A°; e) = V^y.e -f X° . (7.4) 

where Y°-'',u'j^ denote the collective coordinates and recoil velocities of the group of N D-particles, with the indices 
a,b = 1, . . . N. The reader should also note that, as explained in [1^, the range of the variable ui € (0, oo), as a result 
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of the impulse approximation of recoil, which dictates that the recoil vertex operators are non trivial only for target 
times after some definite time (taken to be zero for concreteness in the example discussed here p.Tp V The operator 
insertion {Lq — 1)^^ in (|7.ip represents the string propagator on the thin strip of the pinched annulus. 



D 




(a) 

D 




D D 



(b) 

FIG. 3: (a) World-sheet annulus diagram for the leading quantum correction to the propagation of a string state V in a D-brane 
background, and (b) the pinched annulus configuration which is the dominant divergent contribution to the quantum recoil. 

Inserting a complete set of intermediate string states f/, we can rewrite (|7.ip as an integral over the variable 
q = e'^'^^'^ , where r is the complex modular parameter characterizing the worldsheet strip. The string propagator over 
the strip then reads 

A,(z, z') =J2 f 'i^''' l^^(^) ® (ghosts) ® £i{z')\ ^ (7.5) 

where A/ are the conformal dimensions of the states Ej. The sum in (|7.5p is over all states which propagate along 
the long thin strip connecting the discs S and S' (in the degenerating annulus handle case of interest here, S' = E). 
As indicated in (|7.5p . the sum over states must include ghosts, whose central charge cancels that of the worldsheet 
matter fields in any critical string model. 

In ()7.5p we have assumed that the Virasoro operator Lq can be diagonalized in the basis of string states with 
eigenvalues their conformal dimensions A/, i.e. 

Lo\£i) ^ Aj\£j) , q^'>-^\£j) ^ q^'-^\£j) (7.6) 

However, this simple diagonalization fails jl^ in the presence of the logarithmic pair of operators (|7.2p . due to the 
non-trivial mixing between C and D in the Jordan cell of Lq. Generally, states with A/ — may lead to extra 
logarithmic divergences in (|7.5p . because such states make contributions to the integral of the form J dq/q ~ log 5, 
in the limit q S —^ representing a long thin strip of thickness S. We assume that such states are discrete in the 
space of all string states, i.e. that they are separated from other states by a gap. In that case, there are factorizable 
logarithmic divergences in (|7.5p which depend on the background surfaces S and S'. These are precisely the states 
corresponding to the logarithmic recoil operators (|3.9p and (|3.16p . with vanishing conformal dimension = — 
as e ^ 0+. 

In the case of mixed logarithmic states, the pinched world-sheet topologies of fig. [T] are characterized by divergences 
of a double logarithmic type which arise from the form of the string propagator in (|7.ip in the presence of generic 
logarithmic operators C and D, 

Jdqq^^-^{C,D\(l ^°f'^)|C,Z?) (7.7) 
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As shown in '2(J\, the mixing between C and D states along degenerate handles leads formally to divergent string 
propagators in physical amplitudes, whose integrations have leading divergences of the form 



dq 



log 9 / d^z D{z-e) / d^z' C[z';e) ~ {logSf / d^z D{z;e) / d^z' C{z';e) 



These (log (5)^ divergences can be cancelled by imposing momentum conservation in the scattering process of the light 
string states off the D-brane background [11] . 

We notice at this stage that for brevity we have discussed in this work the simplistic case of a single D-particle 
interacting with an (open) string state. However, as a result of flux conservation, isolated D-particles do not exist, as 
we have mentioned above, and hence in realistic situations one deals with groups of N (with N varying) D-particles, 
interacting among themselves with the exchange of stretched flux-carrying strings. Such groups of N particles can 
be represented by a non-Abelian Wilson loop operator for the gauge group SU (N) . Within the framework of the 
auxiliary field representation of the Wilson loop operator [l^ , the effective abelianization of the matrix cr-model leads 
to a relatively straightforward generalization of the proof of the cancelation of modular infinities by introducing recoil, 
as we now demonstrate. 

The pertinent bilocal term induced by (j7.8p . which exponentiates upon summing over pinched topologies, can be 
written as a local worldsheet effective action using the wormhole parameters [pc,D]i^ to give [isj 



lim 



dpc dpD exp 
WsIpc] 



N 

E 

a, 6=1 
ab pi 



N 



2g2(log<5)- 



+ 



27ra' 

igs [pd] 



'0 

lb /-i 



ds C{X'>{s);e)Us-£)Us):r'''{s) 

ds 



2na' 



ds D{X"{s);e)Us-e)Us) T^'is) 

ds 



(7.9) 



Here we have for simplicity considered only the zero frequency modes of the fields involved with respect to the 
Fourier transformations defined in [ll|. They will be sufficient to describe the relevant cancelations. In (|7.9p the 
(dimensionless) moduli space metric G^*^ (where L, M = C, D) is an appropriate off-diagonal 2x2 matrix, which is 
required to reproduce the initial bilocal operator with the CZ?-mixing of the logarithmic operators. This off-diagonal 
metric includes all the appropriate normalization factors A/l for the zero mode states. These factors are essentially 
the inverse of the CD two-point function, which is finite. 

We consider the propagation of two (closed string) matter tachyon states T12 — e^'^'^^'^-''^' in the background of 
(|7.9p at tree level. In what follows the effects of the C operator are sub-leading and can be ignored. Then, we are 
interested in the amplitude 



N 



A. 



CD 



ec'(l) 




ic' = l 

= ^lim J2 J '^Pc dpD j Dx D^ D^ I 

c' — l 

xexp -N'^So[x]-Y^ 



'(0) 



X Ti [x] T2 [x] exp 

, ^9s[pD\t' 



c=l • 
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a,b=l 



f ds S,c{s - e)—^c{s) 



■G"^" E G%APL\l'[PM]f 



c,d=l 



27ra' 



ds D{X°{s); s) U{s - £)6(s) ^x'{s) 




Cc'(l) 



(7.10) 



where . . . represent sub-leading terms. The scaling property (|3.27[) of the logarithmic operators must be taken into 
account. Under a scale transformation on the worldsheet the C operator emerges from D due to mixing with a 
scale-dependent coefficient a/o' t. This will contribute to the scaling infinities we are considering here. 

The composite D operator insertion in (|7.10p needs to be normal-ordered on the disc. Normal ordering in the present 
case amounts to subtracting scaling infinities originating from divergent contributions of D{X'^{s);e) as £ ^ 0"*". To 
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determine these infinities, we first note that the one-point function of the composite D operators, computed with 
respect to the free cr-model and auxiliary field actions, can be written as 



\c' — l \a,b— 1 

^^-.,(0)expU Y: n^' / {D[X\s);e)D[X\s'Ye))^ 

c' = l \ a,b,c,d { T^Ct ) Jq 



(7.11) 



where we have used Wick's theorem. The second equality in (j7.1ip follows after removing ambiguous 6(e) type terms 
from the Wick expansion in the auxiliary fields using the renormalization scheme described in appendix B of p^ . 
One finds that this procedure has the overall effect of replacing the product of auxiliary fields in the first equality in 
(|7.1ip by the delta- functions SadSbc- 

In what follows we shall ignore, for simplicity, the basic divergences that come from the fundamental string propa- 
gator in (|7.1ip . Such divergences will appear globally in all correlators below and will not affect the final result. As 
a consequence of the logarithmic algebra and the scale transformation p.27p , there are leading (scaling) divergences 
in (17. lip for e ^ O''" which behave as 

g'M-'^'ttr[pnUPDY (7.12) 

Thus, normal ordering of the D operator amounts to adding a term of opposite sign to (|7.12p into the argument of 
the exponential in (|7.10p in order to cancel such divergences. 

Let us now introduce a complete set of states into the two-point function of string matter fields on the disc, 

(TiT2)o = J2 l-^^l' (^il^^)o (^^1^2)0 (7.13) 

where Afj is a normalization factor for the fundamental string states (determined by the Zamolodchikov metric). 
Taking into account the effects of the C operator included in D under the finite-size scaling (I3.24p . we see that the 
leading divergent contributions to (|7.13p are of the form 

(TiT2)o ^ ~V^t (Ti|C)o {C\T2)o + ■■■ (7.14) 

where we have used (13.201) . as well as the fact (c.f. p.l7p . p.lSp and p.l9p ) that the Zamolodchikov metric in the 
C, D basis behaves as [2Q| : 

Gcc ^ , Gdd ^ , GcD ^ Gdc ^ const. . (7.15) 

We now notice that the C deformation vertex operator plays the role of the Goldstone mode for the translation 
symmetry of the fundamental string coordinates x*, and as such we can apply the corresponding Ward identity in 
the matrix cr-model path integral to represent the action of the C deformation on physical states by —iS/Sx"^ [li^. i2(]|. 
The leading contribution to (|7.13p can thus be exponentiated to yield 

^ f ( ^ d \ 

{TiT2)o ^ ^hm E y ^'='(0) exp f -A^'^oM ~J2J^ ds Us - £)^^c(s)j 

' ds ds' Us - e)a(s)6(s' - e)Us') 




6xi{s) Sx'^{s') 



T2[x] ec'(l) (7.16) 
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where we have used the on-shell condition Tj{-^-^)Tk = for the tachyon fields. (|7.16p expresses the non-abelian 
version of the Ward identity in the presence of logarithmic deformations. 

Using ()7.12p , ()7.16|) and normalizing the parameters of the logarithmic conformal algebra appropriately, it follows 
that (|7.10p can be written as 



AcD = ^lim I] jdpcdpD J Dx ^,,{0) 

c' — l 

X exp (^-N^SoM ~ E / ^"-(^ - ^)^^-(^)) 

" ^ f 1 LA/ ^ r b 

^ 1 ~2a^(los6)^ ^ GZ,cd[pL]f[PM] 

a, 6=1 [ WsV B ; ^^^^^ 
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6xi{s) 
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5s Jo 



ds 5(s - S)^a{s) 



dxj{s) 



T2[x] ec'(i) 



(7.17) 



From (|7.17p it follows that the limit t ^ oo localizes the worldsheet wormhole parameter integrations with delta- 
function support 



N 9 



,.fc=l 1=1 



(7.18) 



where (fci,2)i are the momenta of the closed string matter states. This result shows that the leading modular di- 
vergences in the genus expansion are cancelled by the scattering of (closed) string states off the matrix D-brane 
background. Upon rescaling by g^, averaging over the auxiliary boundary fields, and incorporating (|7.18p as an 
effective shift in the velocity recoil operator (c.f. 



(7.19) 



viewed as position operator in a co-moving target space frame, we can identify this renormalization as fixing the 
velocity matrix 



ab 



ab 



(7.20) 



of the fat brane background. Thus momentum conservation for the D-brane dynamics guarantees conformal invariance 
of the matrix u-model as far as leading divergences are concerned. 
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Appendix B: Non-Extensive (Tsallis-type) statistics 



In this Appendix we review the analysis of [27| in constructing classes of stochastic differential equations with 
fluctuating frictional forces which generate dynamics described by non-extensive statistics in the sense of Tsallis [lB| ■ 
Our starting point is the linear Langevin equation describing the dynamics of a Brownian particle: 

u = ~-fu + aL{t) (7.21) 

where L{t) is a Gaussian white noise, 7 > is a friction constant coefficient, and the (real) parameter a described the 
strength of the noise. 

We can take the probability density of the velocity field u to be Gaussian with average {u) = and variance 
(w^) = with /3 = ^ , which plays the role of an inverse temperature of the Brownian particle, taken for simplicity 
to be of unit mass. To get Tsallis statistics, we have to allow for the parameter /3 to fluctuate, which can be achieved 
by allowing in whole generality the parameters 7 and /? to fluctuate as well. 

The case studied in 27] . and used in our work in this article, is the one in which the parameter f3 is distributed 
with degree n, i.e. the probability density of /3 is given by: 

" \ mi/2-l__.„ / nf3 



with Po a constant, playing the role of the average of the /3 fluctuations. 

There are many examples which lead to the distribution of the form (|7.22p . for example [13] if P is given by a sum 
of squares of independent Gaussian random variables Xi ,i = 1, . . . n: 

n 2 

/? = ^ , with ip) = Po , iP^) ~ Pi = -Pi . (7.23) 

1=1 

If the time scale over which P fluctuates is much larger than the typical time scale of order 7^^ that the Langevin 
system (I7.2ip need to reach equilibrium, then the conditional probability p{u\P) that the velocity takes on the value 
u, given a value of P) is Gaussian and reads approximately: 



piu\P)^^^expi^--pu'j . (7.24) 

The probability to observe both a certain value of u and P reads then: 

piu,p)^p{u\p)p{p) (7.25) 

and therefore the probability of observing a value u regardless of the value of /3 is: 

p{u) = J p{u\P)p{P)dp , (7.26) 

which can be evaluated (2^, using (|7.22p . to yield: 

Expression (j7.27|) coincides (up to irrelevant proportionality constants) with the generalised canonical distribution of 
non-extensive statistical mechanics of Tsallis [15|] : 

« : ^ . , (7.28) 



(l + |(g-l)u2) 



provided one makes the identifications q = \ + -^^^ , and P = -^zr^Po 

In [27] more general cases have been considered, where the Langevin equation (|7.21|) is extended to include arbitrary 

du ^ 



frictional forces of the form F{u) = —-S-V{u) 
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